
ANALYSIS OF FREE-MOLECULE FLOW 
WITH SURFACE DIFFUSION 
THROUGH CYLINDRICAL TUBES 


by Thaine W. Reynolds and Edward A. Richley 


Lewis Research Center 
Cleveland, Ohio 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 



TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 



□ DT'lBSb 


ANALYSIS OF FREE-MOLECULE FLOW WITH SURFACE 

DIFFUSION THROUGH CYLINDRICAL TUBES 

By Thaine W. Reynolds and Edward A. Richley 

Lewis Research Center 
Cleveland, Ohio 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — Price $2.00 



ANALYSIS OF FREE -MOLECULE FLOW WITH SURFACE 


DIFFUSION THROUGH CYLINDRICAL TUBES 
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SUMMARY 

The integrodifferential equation governing flow through cylindrical tubes 
under combined free-molecule gas flow and surface diffusion has been solved 
numerically for tube length- to -radius ratios of l/l6, l/8, 1, 4, and 16 and for 
a wide range of flow parameters. Several transmission probability relations are 
developed that permit various comparisons of flow with and without surface dif- 
fusion for given tube length- to -radius ratios. General results are tabulated 
for inlet and outlet surface coverage fractions , gradients in surface coverage 
fraction, and transmission factors. 

Specific solutions for tube flows matched to inlet and outlet surface 
flows are given. Results are presented for the single tube solutions only, but 
relations are given by which these results can be applied to models of porous 
media as well. In the matched solutions, the total transmission factor in- 
creased with increasing values of the surface diffusion parameter. The solu- 
tions are shown to be asymptotic to the long tube approximations. 

Results are compared with those of other investigators and differences are 
found to be primarily due to variations in basic assumptions with respect to 
the model. 


INTRODUCTION 


Flow of gases at low pressures through tubes of small diameter (of the 
order of microns) may be influenced considerably by the process of surface dif- 
fusion of adsorbed molecules along the walls of the passage. These surface 
diffusion effects have been of concern in the evaluation of the flow through 
small single tubes (or orifices) such as utilized in the Knudsen cell measure- 
ments of vapor pressure (ref. l), as well as in the study of flow through vari- 
ous porous media. Examples of the latter include powder beds for gas adsorp- 
tion (ref. 2), porous plugs for cathodes (ref. 3), and certain ion rockets 
(refs. 4 and 5). The theoretical approach to treatment of these flows is 
usually a single tube model. 



Free-molecule flow of gases (i.e., flow at very low pressures) through 
passages of various cross sections appears to be adequately treated (refs. 3, 

6, 7, and 8); however, the combined mechanism of free-molecule and surface dif- 
fusion flow has received much less attention. Usually, efforts to solve the 
combined flow problem are directed toward specific applications (e.g., refs. 1, 
5, 9, and 10). The resulting solutions are often restricted to a particular 
system (i.e., to a specific configuration or to a specific gas-metal combina- 
tion) because of the boundary conditions and simplifying assumptions that were 
made. For example, the orifice flow treatment in reference 1 is restricted to 
length-to-radius ratio tubes much less than unity. Treatments such as de- 
scribed in references 9 and 10 are restricted to long tubes (length-to-radius 
ratio >50) because of the use of the limiting form of the Clausing factor and 
other restricting assumptions. The analysis presented in this report, is not 
restricted in length-to-radius ratio, but covers the range from zero to lengths 
where the long tube approximations become applicable. 

The problem is such that more than one equation may be required to de- 
scribe the flow; for example, one for flow through the tube, one for flow along 
the upstream surface, and one for flow along the downstream surface (see 
sketch (a)). The equation describing the combined free-molecule and surface- 
diffusion flow through the tube is the most complex. It is a second- order 
integrodifferential equation not amenable to closed-form solution. Various 
boundary conditions may be specified, which in turn must be matched to condi- 
tions at the inlet and exit ends of the tube to obtain a complete solution. 
These boundary conditions will be discussed in detail in terms of the physical 
requirements of the problem. Transmission coefficients for the combined flow 
are also developed and presented. 

An important feature of the tube flow solutions presented herein is that 
they are obtained by an iterative numerical method that exhibits very rapid 
convergence and is completely general. The technique is described in appen- 
dix D by Carl D. Bogart. 


ANALYSIS 
The Model 

Establishing an appropriate mathematical model from the physical one ap- 
pears extremely difficult if all possible aspects are to be included, that is, 
if it is to represent a pore in a piece of porous material as well as a simple 
cylindrical hole in a plate. In the former case, where gross flow characteris- 
tics are primarily of interest, recourse is often made to the inclusion of un- 
certainty factors applied to macroscopic flow relations without regard for the 
detailed physical model (ref. 2). These factors must then be evaluated experi- 
mentally for specific systems. Another macroscopic approach (described in de- 
tail in ref. 4) treats the porous medium as if it consisted of parallel capil- 
laries having some mean capillary radius. Physically this model also seems to 
be satisfactory for a microscopic approach as long as it is recognized as merely 
a convenient tool for mathematical analysis. The analysis can be constructed 
around a single typical capillary, and it is equally applicable to a wide variety 
of problems. The model used herein is depicted in sketch (a), which shows the 


2 




gas -phase and surface -diffusion compo- 
nents of the flow. Throughout this 
report the use of the term gas-phase 
is restricted to mean that component 
of particle flow in free space and 
does not include any particle flow 
along surfaces. That flow which in- 
cludes both components will always 
have the term total attached to it. 

Ahead of the upstream end of the 
tube, the gas is in equilibrium. Par- 
ticles arrive at and evaporate from 
the wall leaving it with a net surface 
coverage. Because of surface diffu- 
sion, the surface coverage is shown to 
decrease along the upstream surface in 
the direction of the tube entrance. 
Surface-diffusion flow takes place 
through the tube also, and in this 
case there is an accompanying gas- 
phase flow. The gas -phase flow is 
assumed to be free-molecule flow; that 
is, particle -particle collisions are 
negligible and only particle -wall col- 
lisions must be considered. Along the 
downstream surface, which is exposed 
to vacuum, flow is by surface diffu- 
sion with evaporation. Thus, it can 

be seen that the system can be broken into three segments. These segments will 
be described in detail in the next section and the governing mathematical rela- 
tions will be developed. Additional appropriate assumptions or restrictions 
regarding the flow will be introduced as they become applicable. 
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Cross section through single capillary 
(a) 


The Tube 

Flow through the tube consists of both gas -phase and surface-diffusion 
flow. Particles arrive at and depart from the tube wall in much the same way 
that takes place on the upstream surface; however, as will be seen, the mathe- 
matical relation describing the process differs considerably. Particle leaving 
rate from the tube wall when no surface-diffusion effect is considered is given 
by the following equation first derived by Clausing (ref. 6) and developed in 
detail in reference 8: 


v(x 2 ) = 2 + f v(x)F 3 2 dx (1) 

’ 0 

(Symbols are defined in appendix A.) The barred quantities are dimensionless 
ratios with respect to the tube radius R. Sketch (b) will aid in visualizing 
the physical significance of the terms in equation (l). 
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Location 1 


2 


3 



The left side of equation (l) is the local 
net leaving rate at an element of area at lo- 
cation 2 and is assumed equal to the local net 
arrival rate. On the right side of the equa- 
tion, the term qqFi 2 ■ is arrival rate of 
particles directly from the inlet. The quan- 
tity Fj_ g is a geometric factor arising from 
the assumption that the flow has a random di- 
rectional distribution (ref. 8): 


■ 1,2 


5^ + 2 


+ 4 > 


1/2 


- 


( 2 ) 


The second term on the right side of equation (l) is the arrival rate of parti- 
cles from the tube wall (location 3) integrated over the complete area of the 
wall. It will be recognized that v(x) is the function of which v(x£) is a 
local value. The geometric factor in this term is 


P 


3,2 




(x - x g ) 2 +6 

[a - + 4 f /2 


• ► 
J 


(3) 


An analytic solution to equation (l) has not yet been published; however, 
numerical solutions have been obtained by several investigators (refs. 8, 11 
to 19). The results show that v(x) is almost linear. An interesting feature 
of this near linear relation is discussed in appendix B. 


If an adsorbed layer is present on the tube walls, then a term accounting 
for surface-diffusion flow must be added to equation (l): 


v(x 2 ) 


“ ^ 1 F 1,2 + 


/ L 

v ^) F 3,2 ^ + °"( x 2 ) 


( 4 ) 


This added term represents the net accumulation of particles at location 2 due 
to flow by surface diffusion. It is the solution of equation (4) that is of 
primary interest to describe the flow through the tube. (A derivation of 
eq. (4) is given in ref. 3.) Before a solution of equation (4) can be at- 
tempted, a relation between the local leaving rate v and surface concentra- 
tion a is required. Where experimental data are available for a specific 
system they may be used (see ref. 5). More generally, various analytical rela- 
tions for surfaces at constant temperature (adsorption isotherms) have been 
developed and are given in the literature. A detailed discussion of all these 
is beyond the scope of this report; however, appendix C is included to give an 
insight into some of the fundamental assumptions and common relations. As dis- 
cussed in appendix C, the relation used herein, which is also an appropriate 
approximation of all the relations given, for low values of surface coverage 
fraction Q(Q = °/v m )> is 
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V 




= c^e 


(5) 


A useful concept is the time required to deposit a monolayer onto a clean 
surface assuming no evaporation occurs: 


T 


m 



( 6 ) 


Then, when the evaporation rate v equals the arrival rate \jl the surface 
coverage fraction can be designated (from eqs . (5) and (6)) as 


6 


CO 



T 



(7) 


The surface coverage fraction 0 M is the maximum value of coverage that is 
expected in any system in the sense that it is at the maximum arrival rate and 
occurs in the absence of surface diffusion effects. Then, at any other evapo- 
ration rate v(x), again from equations (5) and (6), 

a 

v(x) a % Q m _ 0(x) (g\ 

LU T a T 0 v ’ 

K 1 m — oo 

T 

m 

Thus, by using equation (8), equation (4) can be expressed in terms of 0 as 

0 ( x 2 ) = F l,2 + r 0(x)P 3 g dx + CgS'^Xg) (9) 


where 


e(x) 


e(x) 


do) 


and 


C 


2 




m 




2 6 oo 


( 11 ) 


With the specification of Cg, L and two boundary conditions, equation (9) 
can be solved by using iterative numerical techniques. The required boundary 
conditions £an_be specified by selecting any two of the variables 0(0), 0(L), 

0 ' (0), or 0*(L). The choice is a matter of convenience with respect to the 
particular numerical procedure used. In the one used herein the values of 
0(0) and 0(L) are specified as boundary conditions (see appendix D). Values 
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of Q x (needed later) are then obtain- 
able from the resulting solutions. 

Mathematical solutions of equa- 
tion (9) may be obtained without regard 
for upstream surface or downstream sur- 
face flows. As will be seen however , 
depending on the physical model of the 
system, only certain of the solutions 
may be applicable . 


The Upstream and Downstream Surfaces 

For the model chosen as illustrative of one physical system the upstream 
and downstream surfaces may be regarded as annular surfaces of inner radius R, 
as shown in sketch (c). 

Upstream surface . - The differential equation describing the steady-state 
surface concentration a is (ref. 5) 

^| o "( r ) + - i - cr ' ( r )1 = R 2 (v - jjl ^) ( 12 ) 

L r J 

Equation (12) can be rewritten by using equations (5) and (ll) as 

r 2 e"{ r) + r0'(r) - r^e = -C^r 2 (13) 


where 


C 


_ R^ _ _1_ 
3 ~ ~ C 2 


from equations (6), (7), and (ll). 


(14) 


Equation (13) is in the form of a modified Bessel equation for which the 
general solution is 

e{x) = bi 0 (?Vci) + cK 0 (?VC3) + 1 (15) 

For a single tube or hole in a plate, as r gets large, 0(r) approaches 0 ro . 
This requires the constant b _to be zero since Iq goes to infinity and Kq 
goes to zero with increasing r. For this case 


e(r) = cK 0 (r-v/cJ) + 1 


(16) 


In general, however, for closely spaced tubes (pores) the constant b is 
not zero. Assume that at some maximum radial distance p the concentration 
gradient is zero; that is. 
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e'(p) = 0 


(17) 


Applying this boundary condition to equation (15) yields 

b _ %(P 

c " ixCp-Vci) 

By using this relation, then, equation (15) becomes for the general case 

Kn(pVC^) r—s 

! - &38f> Io(rV5?) * K ° (rVC5> 


(18) 


1 - 0 T 


Ki(p-v/C^) 

Ii(pV%) 


(19) 


I o(V^3) + k o(V c 3) 


For widely spaced pores this reduces to the single tube relation (i.e., for 
P -* “ ) ) 

1 - 0(r) _ K p( r V^3) 

1 - 0 R k 0 ( V^i") 

Differentiation of equation (19) yields for the slope 


(20) 


0 ! (r) - - ( 1 - 0^ ) V ^3 


K-i ( p -yj C, ) 

I x (pV%) I l (r ^3 - K l( r Vci) 


Kp( p -y/C 3) 


( 21 ) 


Ii(pV^) 

which for the single tube at r = R becomes 

e R = (1 " 


Iq(VS ) + K o(VcJ) 


Ki( -v/Cj) 


VVc? 


( 22 ) 


With Cg specified equation (22) may be solved for various values of 0 R . 
This equation will be considered in detail in the section Matched Boundary 
Solutions . 

The quantities 0 M and C 3 are not independent but are related from 
equations (7) and (14): 


C 3 = K 2 ^^)- 1 


(23) 


Downstream surface . - There is no gas-phase arrival rate on the downstream 
surface (see sketch (c7), and the differential equation describing the steady- 
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state surface concentration on this surface is 


9 


Jjr"(r) + -i- a'(rjT| = v^R^ 


( 24 ) 


A development similar to that used for the upstream surface leads to the 
following solutions : 

For the general case of closely spaced pores 

k i(pV%) — - — — , — . 

^ I 0 (rV%) + 


0(r) _ 

e R K^pVcj) 

and for the single tube or pore 


io(vrc£) + k 0 (Vc^) 


9(r) = K o( r V^i) 

% k 0 (vc;) 


The slope for the general case is 

K i(pycj) 

i^pVcJ) 


Il( r V^3) - Ki(r-\/C3 ) 


0'(r) = 0 R VCJ 


Ki(pVci) 

^(pVcJ) 

and for the single tube at r = R 


io(Vc^) + Ko(Vci) 


Ki(Vci) 


Kq( V%) 


(25) 


(26) 


(27) 


(28) 


With C 3 specified equation (28) may be solved for various values of 0 r. 

This equation will be considered in detail in the section Matched Boundary 
Solutions . 

The similarity between the upstream and downstream relations is apparent. 
Figure 1 is a plot of the specified functions that appear in both equations (21) 
and (27). This combination of functions is herein defined as 
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jr(c 2 ) = 


Ki(FVcj) 

j^pysp 

Kg(pVjg) 

I-l(pVS) 


l!(rVC^) - Ki(rVc^) 


I 0 (Vci) + Ko(Vci) 


(29) 


and is plotted against the variable C£ in figure 1, 0 % being the reciprocal 
of C 3 (eq. (14)). 

It is apparent from figure 1 that for values of Cg less than about unity, 
the 0* solutions for the single tube case (p = oo) are applicable to closely 
spaced pores (p ** 2) as well. Subsequent solutions will be presented only for 
the case p = °o. The solutions for cases of p less than infinity may be ob- 
tained by following the same procedures to be outlined for the single tube case, 
but by using the appropriate curve from figure 1. 


fetching Boundary Conditions 

Three equations of primary interest for the physical model are equation (9) 
describing the tube flow, equation (2l) describing the upstream surface flow, 
and equation (27) describing the downstream surface flow. These equations are 
mathematically independent in that each may be solved with no knowledge of the 
others. The resulting solutions may be combined to represent a physical system 
provided some interrelation is specified. 

In the treatment of the physical model given herein, complete solutions 
that satisfy the following boundary conditions will be presented: 


8q = 8p> (upstream) 

(30a) 

0* = 0* (upstream) 
0 R 

(30b) 

0 L = 0 R (^°™ s bream) 

(30c) 

0 =0 (downstream) 

L R 

( 30d) 


In addition, a wide variety of independent tube flow solutions will be given 
from which it is possible to construct complete solutions to other physical 
systems . 

The same boundary conditions (eqs. (30a) to (30d)) are specified in refer- 
ence 5. In reference 1, 0 r at the upstream end is equated to 8^, thus elim- 
inating all upstream surface diffusion effects from the system. At the down- 
stream end the required matching condition of reference 1 is the same as equa- 
tion ( 30d ) • 
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Transmission Factor 


The various transmission factors that are 
defined in the case of gas -phase transport 
alone are all ratios of the net flow out the 
exit to the inlet gas -phase arrival rate. For 
that case there is no other flow mechanism 
introducing material into the tube; however, 
with surface transport possible, as in the 
present case, some mass flow may be introduced 
into the tube by surface transport also. It 
may, in fact, exceed that introduced by the gas-phase transport. Nevertheless, 
it is still convenient to define transmission factors based on the inlet gas- 
phase arrival rate alone, even though the transmission factor may then exceed 
unity. The inlet gas arrival rate is a readily evaluated number, independent 
of the tube flow conditions, while the surface transport portion of the inlet 
flow is not. The various transmission factors presented herein, therefore, are 
all referred to the inlet gas arrival rate only. 


Surface- 
diffusion 
flow, n 2 
Gas-phase 
flow, nj 


V 

\y////////////ZA 


J. 


n 4 

n 3 


^777777777777^3 

/ A 

Inlet Exit 

plane, plane, 

x * 0 x * L 

<d) 


Consider the system depicted in sketch (d). The mass flow (particles/sec) 
crossing the inlet plane due to gas-phase flow is 

n l = ^i^ 2 (31) 

and that due to surface-diffusion flow is 


n 0 = -2k 




dx/. 




(32) 


At the exit plane the mass flow expressions are more complex. The gas-phase 
flow is 


n^ = Pgi-L^kR 2 (33) 

where Pq is a gas-phase transmission factor, defined as the ratio of the gas- 
phase flow passing the outlet plane to the gas -phase flow entering across the 
inlet plane. As shown in reference 8, Pq. may be considered to consist of two 
components, a direct and an indirect transmission factor: 

P G - p a + p i < 34 > 

where P^ is simply the fraction of entering particles that pass directly 
through the tube without any encounters with the wall. It is a function only 
of L: 


P 


d 



(35) 


The flow Pj_ is the ratio of the exiting gas-phase flow coming off the tube 
wall to the entering gas -phase flow. It is given by the following integral 
(ref. 8): 
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e(x) 


'0 ~0 

The function 0(x) is given by equation (9). 

The surface -diffusion flow out the exit is 


(L - x)Ql - x) 2 + 1 - 
|[(L - x) 2 + 1 + r 2 ] 2 - 4r 2 }- 5 ^ 2 


dx dr 


n 4 


■ 2irR ^ (H) L = - 2lt ^ a m 0 L 


The total transmission factor is, then., 

n 3 + n 4 


P T = 


n l 


= P G + P S 


where Pg is a newly defined surface transmission factor 

n 4 

P Q — = -2C_e T 

S nj_ 2 L 


(36) 


(37) 


(38) 


(39) 


Other Considerations 


General numerical solutions of the tube_ flow equation (eq. (9)) are ob- 
tainable with the specification of 0%, L, 9 0 , and 0 L . Since solutions are 
obtained by an iterative procedure, specification of a convergence factor e 
is also required. Convergence is defined herein in terms of this factor as 


0 ( x )n+l ~ e t*)n 

e(x) n 


< e = 10“5 


max 


(40) 


Before going on to the general and matched-boundary solutions, a brief 
discussion on the makeup of the 6 m and C 2 terms may aid in considering the 
physical interpretation of the results. As seen from equation (7), specifying 
9 from physical considerations requires a knowledge of the upstream equilib- 
rium arrival rate pq and the parameter Cq. For a given gas pq is related 
to the pressure and temperature : 

P x = p(2jtMkT) -1 / 2 (41) 

As an example, figure 2 is included to show the variation of pq with T for 
a few substances, where the vapor pressure relations of p with T from 
reference 20 have been used. The parameter Cq is expressed in terms of the 
quantities a , r, t q , Q, and T gV from equation (CIO). Some of these are 
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fairly constant, but others are highly uncertain and may vary widely with mate- 
rial. For many materials the monolayer coverage o m is on the order of 10 15 
atoms per centimeter squared. As mentioned in appendix C, the sticking coeffi- 
cient x has been assumed to be unity. The quantity Tq is a proportionality 
constant relating the mean adsorption lifetime T with the desorption energy 
Q and the temperature of the system (eq. (C7)). A reasonable value of t q is 
about 1C>~13 second although values considerably larger than this have been 
reported. Based on the value of 10" figure 3 gives the variation of log^QT 
with Q for various values of T. If the relation between Q and T is 
known for a given system, all the quantities necessary for evaluating 6 are 
available . 

As shown in equation (ll), two additional parameters are required to eval- 
uate Cg: the surface- diffusion coefficient 3 and the tube radius R. Ex- 

perimental values of 3 are very sparse. A value of 10"^ centimeter squared 
per second was assumed for silver vapor on molybdenum and on nickel in refer- 
ence 1. At high temperature the value of 3 for cesium on tungsten is also 
about 10"^ centimeter squared per second (ref. 2l). A theoretical approach to 
describing 3 is also given in reference 21, where 3 is considered to depend 
upon the adsorption time t and the mean transit time t. The relation 
employed Is 


^ 1 a 

® ~ 4 T + t 


(42) 


where the mean transit time across an adsorption site t is 

a a 

t = - = — f7p ( 43 ) 

frkT Y 2 

\ 2m ) 

At low temperature the adsorption time is much longer than the transit time, 
and 3 exhibits an exponential temperature dependence. At high temperature 
the transit time becomes controlling, and then 3 varies as the square root of 
the temperature . This behavior is portrayed in sketch ( e ) . At the higher tem- 
perature levels then, it is reasonable to expect 3 to be of the same order of 
magnitude for different systems, its variation being only inverse with the 
square root of the particle mass. The use of the term high temperature is 
only relative, of course, the important factor being the relation of the ad- 
sorption and transit times. 

This in brief describes the physical parameters of interest. It should 
also be pointed out that with regard to the upstream and downstream surface 
flow relations, (eqs. (21) and (27)), no additional parameters are required to 
define C3. 

The following sketch recapitulates the material discussed thus far. In 
sketch (f ) the circled numbers refer to applicable equations, and the relation 
of p and Q to T is assumed known. As pointed out previously, specifica- 
tion of certain physical parameters imposes interrelations among the C ! s 
(e.g., eq. (23)). Conversely, the C f s may be specified with no knowledge of 
the physical parameters. 
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RESULTS AND DISCUSSION 
General Tube Solutions 


Solutions to equation (9) were carried out on an IBM 7094 computer by us- 
ing the procedure outline d_in appendix D. Initially, solutions were obtained 
for various pairs of 8q, 8^ values for each combination of L and Cg. It 
is not practical to present all the tube-wall coverage data that were obtained. 
Sketches (g), (h), and (i) are shown as typical examples for various levels of 
coverage at the end points. The trend of the solutions which was obtained in 
all cases can be noted from these sketches. When 8 q and 8-^ are fixed, the 
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0 .2 .4 .6 .8 1. 0 


0L 

(0 



1.0 .8 .6 .4 .2 0 


(k) 

solution varies from a straight line 
through 0 q and 0^ at high C 2 /L 

values; and tends toward the straight 
line through A and 1 - A as C 2 /L 
approaches zero (see appendix B), ex- 
cept for holding to the imposed end- 
point values. 

From various plots of the solu- 
tions; it was found that all values of 

0q; 0^, and Pj_ could be represented 

by "maps" as shown in sketches (j), 

(k) ; and (l). These maps can_be con- 
structed from three sets of Qq, 0 ^ 

conditions. These latter results are 
summarized in tables I and_II. Table I 
has values of the slopes 0 q and ©£ 

for the combinations; ( a.)_ 0q = 1.0; 

0^ = 1.0; (b) 0 L = 1.0; 0 O = 0; and 

(c) 0 q = 0; 0^ = 0. Table II contains 

values of the indirect transmission 
factors Pj_ for the same sets of pa- 
rameters . 
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The three values for each C^/L that are given in the tables can be used 
to map the solutions for any set of end-point conditions as indicated in 
sketches (j), (k), and (2). These 3 points in table I are designated a, P, 
and 5 to assist in showing their use in constructing the solution maps. Note 
that the graphs of sketches (j) and (k) are identical when the coordinates are 
labeled in the manner illustrated. 

The gas-phase transmission factor Pq (eq. (34)) differs from Pj_ only by 
the factor P^ and is not tabulated. The direct- transmission probability P^ 
is only a function of L (eq.. (35)) and is shown plotted in figure 4. Similarly, 
the surface-diffusion transmission factor P s is readily obtainable from the 
tabulated values of slope through equation (39) and thus is not tabulated. 

It will be noted in sketches (g) and (h) that for a given L as Cg de- 
creases, the absolute magnitude of the slopes at the end points increases; 
however, the product C^d x , which is a factor in the surface-diffusion flow 
transmission factor, decreases as 0^ decreases, as would be expected* 

Matched Boundary Solutions 

The matched tube and end-plate solutions are readily obtained by using 
the graphs for the general tube solutions and by superimposing the relations 
for the upstream and downstream surfaces (eqs. (2l) and (27)). The procedure 
is illustrated in sketches (m) and (n). 
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In sketch (m) , the ordinate is the slope at the downstream end of the 
tube^ while in sketch (n) it is the upstream value. From the symmetry of the 
two plots it may be readily shown that the solution must be along the line 
®0 + = so that actually only one map is required to obtain the solu- 

tion. For example, (at r = R) equations (2l) and (27) and the matching rela- 
tions (30) may be combined to give the relation 



The general tube solutions, illustrated by sketches (m) and (n), may also be 
expressed mathematically as 

e L = f l( 0 o) + b l 0 L 

and 

= ^l(-*- " + b l(l “ e o) 

where ^(Sq) and f-j_(l - 0p) are bbe intercepts in (m) and (n), respectively, 
and b]_ is the slope of the lines. Combining these two expressions gives 

Oq i"i(i " + b l^i " ®q) 

d L fl( Q 0 ) + b l e L 

[This expression combines with the first relation to give 

1 - @0 %,) + k]_(l " 

e L fi(0o) + b i e L 

Cross multiplication and simplification of this last expression yield the rela- 
tion 


(1 - e 0 )fi(©o) = " 0 l) 

which further reduces to 


1 - 



since ^(Sq) equals f^_(l - 0 when 0^ = 1 - 0 q. The matched-boundary solu- 
tions thus have the properties 


Q 0 + d L ~ 1,0 
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and 



It should be noted that the di- 
rection in which the coordinates x 
and r have been taken herein requires 
that the matched slopes at the upstream 
end be equal in magnitude but opposite 
in sign. This is illustrated in 
sketch (o) showing a typical complete 
solution. 

The solutions to equation (9), which 
are matched to the end-face solutions 
through the boundary_conditions of equa- 
tions (30) and for p = 00 , are presented 
in table III. A plot of the solutions 
for Gj^ for these matched-flow cases is 
shown in figure 5, which shows 0 ^ var- 
iations with L for various constant 
values of C 2 /L. The corresponding 
values of the surf ace- transmission 


(p) factor P 3 for these same conditions 

are shown in figure 6 . Figure 7 shows 
the variation in the total transmission factor with tube length for various 
constant values of Cg/L. The lower curve of figure 7 represents gas-phase 
flow only. Clearly, the surface -diffusion flow can become a large fraction of 
the total flow as C 2 becomes large. The calculations presented herein are 
also seen to extrapolate very nicely to the long tube limits as calculated 
from the relations of reference 10 . 


The matched-boundary-condition solutions for the single tube (p = 00 ) may 
also be examined with regard to their applicability to porous materials having 
capillarylike pores. As depicted in sketch (o), the slope of the upstream and 
downstream surface solutions approach zero with increasing r. The solutions 
would be expected to be valid for pore spacings that allowed this condition to 
be approached. As was pointed out previously (see the section Upstream and 
Downstream Surfaces), the appropriate curve of figure 1 may be used to obtain 
the exact solution applicable to a given pore spacing. The parameter p would 
be equivalent to one-half the distance between pore centers (sketch (a)). 

The function ^(Cg), plotted in figure 1 , when multiplied by the _ 
yields the slope at R (see eq. (27)). The change in this function when p is 
less than infinity can thus be used to obtain the change in slope, and hence in 
surface flow for closely spaced capillaries. The relation is derived as 
follows with the aid of sketch (p), where the coordinates have been changed 
from the corresponding maps of sketches (m) and (n) by using the matching rela- 
tions (30). 
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The solution 

tions 



is obtained by simultaneous solution of the two equa- 


<«*> 

and 

5 E - 5 + -E 

Equation (44b) corresponds to the relation e 0 + %, 
rameters of sketch (p). 

The solution for the slope is 

•T( c 2 )_ 

( q '\ _ . P _ S 

' R'p ~ ^(C 2 )_ - m ° 

The corresponding solution for p = «> is the same as that for equation 
(45) with ^"(Co)— replacing ^”(Co)_. Thus, the ratio of the slope for 

_ c p=co ° p 

finite p to that for p = » is 

*XC 9 )_ #tc )_ - m 

P C P c p=oo 

= ^(c- 2 )_ jr (c - 2 l_ — ' ( 4I 

p=c0 p 

The values for the slope m are included in table I, and values of the func- 
tion ) are plotted in figure 1. Thus, the solutions for closely spaced 

pores (p less than infinity) can be obtained from the single tube solutions 
(i.e., p" = oo ) and equation (46). 

As an example, the case of C 2 = 10, p = 2, and is presented in the 
following table: 



(44b) 

= 1.0 in terms of the pa- 


(45) 


L 

£Uo) 

(fig. 1) 

m 

(table I) 

(s R ) iT=2 //(e RW 


-D| 

!( 

ro 

p = oo 

(Eq. (44)) 

1 

0.14 

0.68 

-2.015 

0.258 


For this illustrative example , then, the slopes and hence the surface trans- 
mission probabilities (eq. (39)) for the closely spaced pores (p = 2) would be 
smaller than that for the corresponding single tube cases by the factors shown. 
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Comparison of Results with Those of Other Investigators 

In reference 1 two systems, silver on molybdenum and silver on nickel, 
were investigated. The values of the physical parameters cited therein and 
also used herein, are given in the following table (notation of the present 
paper being used in all cases): 


System 

■ 

Q, 

T 0; 

9 , 


eV 

sec 

cm 2 /sec 

Silver -molybdenum 

2.0 

icr 10 

10- 4 

Silver-nickel 

1.5 

io- 8 

io- 4 


The resulting values of Pg, the ratio of surface -diffusion flow at the 
exit to the gas-phase flow at the inlet are compared in the following table: 


System 

T 3 

°K 

L, 

cm 

R, 

cm 

CO 

Pk 

Ref. 1 

Present 

report 

S ilver -molybdenum 

1000 

0.005 

0.05 

0.27 

0.20 

S ilver -nickel 

1000 

.005 

.05 

.15 

.086 

Silver-nickel 

1280 

.005 

.025 

.021 

.017 

Silver-nickel 

1280 

.001 

.025 

.037 

.023 

Silver-nickel 

1280 

.005 

.10 j 

.005 

.0046 


It can be seen that the values of Pg of the present paper are all lower than 
those of reference 1. This is expected since one of the assumptions of refer- 
ence 1 was that 0 q = 0^, which is higher than the value solved for herein, and 
thus leads to somewhat higher flow by surface diffusion in the tube. 

A direct comparison of the results of the present calculations with those 
of reference 5 is difficult to make because of the differences in the forms of 
the adsorption isotherms used. The present results are based on a linear 
adsorption relation (eq. (5)), which, as previously pointed out, is generally 
acceptable in representing most adsorption relations at low coverage. In 
reference 5, the evaporation rate data for the cesium- tungsten system (ref. 17) 

have been used. These data are very nonlinear 
in the range used therein, as can be seen in 
sketch (q). This sketch shows the evaporation 
rate curves for atoms, ions, and for the total 
of ions plus atoms for a surface temperature of 
1400° K. In the absence of an electric field 
to remove ions from the surface, atoms alone 
will evaporate. If a sufficiently strong elec- 
tric field exists at the surface to remove the 
ion current, both ions and atoms evaporate. 
Under these conditions, at low coverage the ion 
component is the major portion of the total 
evaporation rate. To use a linear relation to 
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approximate either the atom or ion evaporation rate, obviously means restrict- 
ing the solutions to cases where 6 m is no greater than about 0.08 for the 
atoms or about 0.01 for the ions. 

The use of the total evaporation rate curves as the adsorption isotherm 
relation in the solution of equation (9), as was done in reference 5, involves 
the simultaneous assumptions that an electric field exists to remove ions but 
that the ion trajectories are not then influenced by this field. 

Thus neither the approach of reference 5 nor the present one appears 
adequate to calculate the case for total ion plus atom evaporation of the 
cesium- tungsten system at coverage fractions above approximately 0.01. 


CONCLUDING REMARKS 

The object of this report was to analyze the flow characteristics of a 
cylindrical tube, or pore, under conditions of free-molecule flow with surface 
diffusion. The integrodifferential equation that describes the flow in the 
tube was solved numerically for a wide range of variables. Matching boundary 
conditions were utilized to obtain solutions that included upstream- and 
downstream- face surface -diffusion flow as well. Transmission factor relations 
were developed that allow quantitative comparisons to be made of the differences 
between flow with and without surface-diffusion effects. 

In general, it was found that with the inlet and exit surface coverage 
fractions held fixed, the shape of the curves describing the surface coverage 
fraction in the tubes varied as follows: with a large value of the surface- 

diffusion flow parameter, solutions were nearly straight lines between end 
points } as the parameter decreased, the curves approached the nearly straight- 
line solutions of zero surface-diffusion flow except for the fixed end points. 

Matched-boundary-condition solutions were found to be obtainable from maps 
of general solutions for a given tube length and surface-diffusion flow param- 
eter for either a single tube or for a closely spaced array of tubes. The 
matched solutions were presented only for the single tube model, but relations 
were derived that permit the application of these results to porous media 
models made up of closely spaced capillaries. Transmission factors were found 
to increase with increased surface-diffusion flow parameter, the increase being 
proportionally much greater with longer tubes. 

Comparison of results with one investigator showed excellent agreement. 

In another comparison, differences were attributable to dissimilar basic 
assumptions or boundary conditions for the model. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, October 21, 1965. 


20 



APPENDIX A 


SYMBOLS 

[Unless otherwise specified, cgs units are used throughout.] 

A parameter defined in eq. (Bl) 

a distance between adsorption sites 

B parameter defined in eq. (Bl) 

b parameter defined in eq.. (15) 

b^ slope of lines of general solution maps (sketches ( l ) and (m)) 

parameter defined in eqs. (Cl) to (C5) 

C 2 parameter defined in eq. (ll) 

parameter defined in eq. (14) 
c parameter defined in eq. (15) 

surface-diffusion coefficient, cm 2 /sec 
F 1 2 geometric factor defined in eq. (2) 

2 geometric factor defined in eq. (5) 

^(Cg) function defined by eq. (29) 
f parameter defined in eq. (C3) 

function giving intercept of lines in sketches (m) and (n) 
g parameter defined in eq. (C5) 

Iq modified Bessel function of the first kind of order zero 

1^ modified Bessel function of the first kind of first order 

K q modified Bessel function of the second kind of order zero 

modified Bessel function of the second kind of first order 
k Boltzmann constant, erg/°K 

L tube length, cm 

M mass of atom or molecule, g 
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m 


n i 

n 2 

n 3 


P. 

i 


P T 

P 

Q 

B 

r 

T 


T 


eV 


t 

v 

x 


ct,P,S 

r 

6 


slope of lines on general solution maps for -which 0 O + ®L = 1*0 (see 
table I and sketch (p)) 

gas-phase flow rate into tube, atoms/sec 

surface-diffusion flow rate into tube, atoms/sec 

gas-phase flow rate out of tube, atoms/sec 

surface-diffusion flow rate out of tube, atoms/sec 

direct transmission factor, fraction of entering gas-phase flow parti- 
cles that exit downstream without experiencing any collisions with 
tube wall 

sum of P^ and P^ 

indirect transmission factor, fraction of entering gas-phase flow par- 
ticles that exit downstream after experiencing one or more collisions 
with tube walls 

ratio of surface-diffusion flow at exit to entering gas-phase flow 

total transmission factor 

pressure, dyne/cm 2 

desorption energy, eV 

tube radius, cm 

radial distance variable 

temperature , °K 

T 

temperature expressed in electron volts, T e y = — j 

time, sec 

velocity, cm/sec 

axial distance variable, cm 

axial distance variable at specific location (see sketch (b)), cm 
designation of general solution points (see table i) 
sticking coefficient 

convergence criterion defined in eg. (38) 
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G 


fraction of monolayer surface coverage 


0 l gradient in fractional surface coverage with respect to nondimens ional dis- 
tance variable 

Q U second derivative of Q with respect to nondimens ional distance variable 


value of surface coverage fraction in tube at inlet 
gradient of surface coverage fraction in tube at inlet 
value of surface coverage fraction in tube at exit 
gradient of surface coverage in tube at exit 

surface coverage fraction on upstream or downstream face at tube radius R 
gradient in surface coverage fraction on either face at tube radius R 
maximum surface coverage fraction attainable in system 
arrival rate, atoms/( cm 2 ) ( sec ) 

arrival rate at conditions existing in inlet chamber, atoms/(cm 2 ) (sec ) 
rate of particles leaving surface, atoms/ (cm 2 ) ( sec ) 
particle leaving rate from downstream face, atoms/(cm 2 )(sec) 
radial distance at which concentration gradient becomes zero 
surface concentration, atoms/cm 2 


a 1 surface concentration gradient with respect to nondimens ional distance 
variable, atoms/cm 2 

a" second derivative of a with respect to nondimens ional distance variable, 
atoms/cm 2 

a m surface concentration for a filled monolayer, atoms/cm 2 
t adsorption time, sec 

time to form monolayer (defined in eq. (6)), sec 
Tq constant in eq.. (C7), equal to 10“ 13 sec in typical calculations 
Superscript: 


Barred distance quantities refer to nondimens ional ratios of the variable 
to the tube radius R (e.g., L = L/R); barred parameters involving Q 
refer to ratios of the variable to the value 0 M (e.g., Gq = Gq/q oq ) 
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APPENDIX B 


APPROXIMATION OF EQUATION (l) 


As mentioned in the report, the solution of equation (l) is nearly linear. 
As a result, it has been approximated by a relation of the following form: 


= A + Bx (Bl) 

^1 


The true solution also has the property 


v(x) v(L - x) 

Hi Pp 


(B2) 


as has been shown to be required based on thermodynamic arguments (refs. 8 
and 11). Using equation (B2), the quantity B can be expressed in terms of A. 


B = 



2A 


Substituting into equation (Bl) yields 


v(x) 

Hi 





(B3) 


(B4) 


The parameter A (the value of v(0)/h 1 ) can be determined approximately 
by assuming equation (B4) to be a solution to equation (l), substituting, per- 
forming the integration, and evaluating the resultant expression at x = 0, 
or x = L. The result is 


A = 


1 



1 + 


g 

zJW+l 

~2 



(B5) 


This result was also noted in references 6 and 19; however, the expression in 
reference 19 differs from equation (B5). The reason for the difference is not 
known. It may be due to a translation error. 

The agreement between the A values obtained from equation (B5) and the 
end-point values obtained from converged numerical solutions (ref. 8) is within 
1 percent as shown in the following table: 
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L 

v(0)/n^ 
(ref. 8) 

A 

(eq. (B5)) 

0.5 

0.6037 

0.6074 . 

1.0 

.6738 

.6736 

2.0 

.7576 

.7574 

4.0 

.8366 

.8369 

8.0 

.8990 

.9010 

16.0 

.9417 

.9450 


The variation of A with L is plotted in figure 8, The values can he used 
with equation (B4), when no greater accuracy is required, to obtain values of 
wall fluxes in the absence of surface-diffusion effects. More accurate values 
for cylindrical tubes as well as for convergent and divergent tubes are avail- 
able in reference 8 . 
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APPENDIX C 


ADSORPTION ISOTHERMS 

The discussion of adsorption isotherm relations included herein is in- 
tended merely to give an insight into the subject. Detailed discussions may be 
found in references 22 and 23. The type of relation utilized for a given 
system is usually governed by the assumptions made "with respect to the follow- 
ing factors : 

(a) Adsorbed film mobility: The classification of an adsorbed film as 

mobile or immobile depends on the relative values of the energy barriers to 
diffusion and to evaporation. If the energy barrier to surface diffusion is 
low compared with the evaporation energy, the particles will migrate for some 
time before they are apt to leave the surface; hence the film is classed as 
mobile . 

(b) Variation of adsorption energy (binding energy) with surface concen- 
tration: Adsorption energy may vary because of surface inhomogeneities. It 

may also vary because of dipole interactions among the adsorbed species. 

(c) Relative surface concentration: Certain forms of the adsorption 

isotherm relation have a filled monolayer as the limiting concentration; others 
do not restrict adsorption to a monolayer but permit multilayer adsorption. 

This distinction also generally depends upon the relative adsorption energies 
of the gas for the substrate and of the gas for itself. 

A summary of some of the various forms for the adsorption isotherm rela- 
tion from reference 23 is as follows: 

Mobile film (no interaction) 


P = Cj9 


and 


P 



Mobile film (interaction) 


p c i i _ e exp (i _ e ~ f0 ) 


Immobile film (no interaction) 


P 


= Ci 



e 


(ci) 


(C2) 


(C3) 


(Cl) 
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I mm obile film (interaction) 


M- = c i y~~q exp(-g0) (C5) 

where p is the arrival rate, C^_ is a proportionality constant, and 9 is 
the surface concentration expressed as a fraction of a monolayer. 

Equations (Cl) and ( Ci) are the simplest of the preceding relations. 
Equation (Cl) is not hounded and permits coverages in excess of a monolayer. 
Equation (C4), usually referred to as the Langmuir adsorption relation (ref. 
23), represents the behavior characteristic of many systems of interest where 
the adsorption energy of the gas for the substrate is considerably greater than 
the adsorption energy of the gas for itself (i.e., the sublimation energy). In 
this case, adsorption is usually restricted to less than a monolayer until gas- 
phase arrival rates (or corresponding gas pressures) approach the saturation 
value corresponding to the temperature of the adsorbing surface. 

Equation (C4) can be derived from a simple dynamic equilibrium approach. 
The evaporation rate is assumed to vary directly with the surface concentration 
o of adsorbed particles and inversely with the average adsorption time t 
(ref. 22): 


v 



The adsorption time t is related to the desorption energy Q: 


(C6) 


T 


t q exp 



(C7) 


The adsorption rate is assumed equal to the product of the arrival rate p, 
the sticking coefficient y, and the area available for adsorption 1-0. 
In equilibrium the adsorption rate equals the vaporization rate so that 


v 


py(l - Q) = 0 J exp 
T 0 



(C8 ) 


or 


where 


p = 




(C9) 


(CIO) 
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From equation (C8) it can be seen that for y = 1 and low values of 6 

v = \i bs CjO (Cll) 

This is identical with equation (5), and, in fact, all the adsorption isotherm 
relations approach this form for low 9. As a result, the solutions of the 
flow relations obtained in this report are based on a relation of the form of 
equation (Cll). They can be considered characteristic of systems having con- 
stant adsorption energy over the range of surface concentration considered and 
for which the isotherm behavior described by equation (Cl) applies without 
limit to 9. For the other isotherm relations they are restricted to systems 
with low 6 values . 
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APPENDIX D 


NUMERICAL SOLUTION OF EQUATION (9) 
by Carl D. Bogart 

Equation (9) can be written in finite difference form as 


e( Xi ) = P^ 2 ( x i) + IT [ F 3,2( x i>°) 0 O + 4F 3,2( x i> h ) 0 l + 

Q 

,, 2 (x i ,2h)e 2 + . . + -§■ [e(x i+1 ) + e(x i _ 1 ) - 2e(x i j] 


2F, 


where the integral has been expanded by Simpson's rule and the second deriva- 
tive by a three point formula. The notation is as follows: Xj_ = ih, where 

i = 0 , 1, . . . , N - 1 ; h = L/(N - l); and F3^ 2 (xi,ih) is the kernel function 
of the integral with e(x^) = 0 j_. This arrangement results in N linear equa- 
tions in N unknowns , which can be written in matrix form as A]_0 = F^g* 

Since the values of 6 at x = 0 and x = L are specified, the integro- 
differential equation has a unique solution. 

By inspection, the matrix Aq_ is tridiagonally dominant because the 
terms on the diagonal, subdiagonal and superdiagonal are on the order of C^/h^. 
The mesh increment h may be chosen small enough so that C^/h^ > 10^ and 
the absolute value of any other terms is less than unity. By the theory of 
regular splittings (ref. 24), the matrix equation may be written as 
(Ti - B-]_)0 = F]_ 2 , where Tj_ is a matrix having only the tridiagonal terms 
and B]_ has no* tridiagonal entries. The iterative scheme then becomes 

T 1 0 n+1 = B-^0 n + F-^ 2 , where n denotes the n^* 1 iterate (n = 0, 1, 2, . . .). 

An initial guess is supplied for 0 and rapid convergence results since B]_0 n 
represents only a small correction and 0 n+ l may be calculated explicity by 
an algorithm for solving tridiagonal systems (ref. 24). 

A schematic flow chart and Fortran IV listing of the computer program are 
as follows : 
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Flow chart 


Calculate 
index of 
0' s to be 
printed out 


M = N - 1 


Read values 
L, D, C, 9n, 
e L- p d 


h * L/100 
Calculate 
multiples of 
h to be used 


Do 10; J g 1, l\T 

/ Calculate ^ 
kernel 

V function i 


D o 90; K = 1, N 

/ Calculate 
entries to 
\ matrix Bi 


Iteration cycle 


Do 136; J =2, N Do 135; J =2, N Do 130; J - l f N 


\ /Intermed. \ / 

Calculate \ / ca l cu .L a ! i f n of \ / Calculate 




'Convergence Not converged / Iteration \ exceeded 
\ test / \ limit test / 


Converged 



Fortran statement numbers 
‘‘‘Fortran do statement 


Do 11; J = 1, N 


Calculate 
F^ 2 and first 
y guess of 0 / 


Do 20; J = 2, M 


Calculate 
diagonal 
of l x 


Do 50; J f J =3, M 


/ Calculate \ 
sub-and super- 
diagonals of 
\ Tj / 















I 


MOLECULffR FLOW PROGRAM 
$ I BFTC REYNQ DEBUG , L I ST , NODECK , R E F 

UIWFNsrON G ( 1 0 1 ) » KER < 101 ) »AA( 101,101 ) »F('l01 ) ,Y( 101 ) ,Al < 101 ) 
DIMENSION WW( 101 ) ,GG< 101 > , A2 < 101 ) ,A3 < 101 ) 

DOUBLE "PRECISION AA,A1,A<:, AB» ‘AC » A. C* DH2* 

2 DH* D, ED, FPS, F. G, _ H6* H, 

3 KER, LBAR, SUM, X J , Y ,EP 

DOUBLE PRECISION GG,WW,SAV,A3 

dimension indexC27> 

DATA N/101/ 

DATA IT, EP/199, 1,0-5/ 

C LBAR=L 
C H=MESH SIZE 

C D = MULTI PLI ER OF D2THET A/DX2 
C D IS ACTUALLY D/C 
C LBAR/H MUST BE EVEN 

C LBAR/H MUST BE LESS THAN OR EQUAL TO 100 
C C IS MULTIPLIER OF y-function 

WR I TE ( 6 » 1 00 ) 

100 FORMAT ( 1 H 1 ) 

C INDEX OF VALUES TO BE PRINTED 
DO 162 J= 1 * 6 
L=2* ( J-l ) 

INDFXI J)=L 

162 ' IITDEXI J+21 ) =L+90 
DO 163 J= 7 ,21 

L = 5* l J-4) 

163 I NDFX ( J ) =1 
M = N — 1 

1 READ (5,105) RLBAR,RD,RC,READ1 ,READ2,PD 

LBAR=RLBAR 
D = RD 
C = RC 

H= • 01 DO*l BAR 
H3=H/3. 

H6_=H/6.D0 
KT 10 J = i,N 
XJ=J-1 
XJ=XJ*H 
AB=XJ**2+A.D0 
AC=DSQRT ( AB ) 

10 KER ( J ) =1 .DO-XJ* ( XJ**2+6.D0 ) / t AB*AC ) . 

C SET UP FINITE DIFF EONS 
DO 90 K=1,N 
JJ=N-K+1 

AA(K»N)=-H6*KER ( J J ) 

A A ( K , 1 ) =-H6*KER ( K ) 

AA ( K , K ) =0 . DO 
L = K— 1 
LL=K+1_ 

1KK.GT.1I AA(K,L)=0.D0 
IF(K.LT.N) AA(K,L)=0.D0 
L = L - 1 
LL=LL+1 

I F ( L • LT ,2 ) GO TO 60 
DO 70 J=2 > L 
XJ = 4-2*M0D< J»2') 

I =K- J+ 1 

70 AA ( K » J ) =-X J*H6*KER ( I ) 

60 IF(LL.GT.M) GO TO 90 

DO "3& "J=LL ,M 
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XJ=4~2*MPD( J»2 ) 

I = J-K+l 

30 AA(K»J)=-XJ*H6*KFR( I ) 

90 CONTINUF 

PH=P/H**2 
PH2=2.P0*DH+1 . D n 
DP 30 J = 3,w 
XJ=4~2*M0D( J,2 ) 

A3 ( J — 1 ) = — DH-X J*H6*KFR ( 2 ) 

XJ=4~2*M0D( J-l »2) 

50 A 2 ( J ) = -DH->XJ*H6#KFP ( 2 ) 

DO 20 J = 2 » M 
XJ=4~2* M PD( J,? ) 

2 n A 1 t J ) =DH2-X J*H6 

A2 ( 2 ) =-DH-H6«^FP ( 2 ) 

A3 t N-l ) =A2 (2 ) 

DO U J=1*N 
XJ=J~1 
X J=XJ*H 
AB = XJ**2 + 4.r>0 
A0= PSQR T { AB ) 

G ( J ) =. ?pp* ( ( XJ**2+2.PO ) /AF-X J ) /C 
11 Ft J ) = ( LBAP+1 .DP-XJ ) / (G*( LP.AR + 2 .P° )) 

C MATRIX FMTRIFS TP ALLOW FNP VALUER TP REMAIN FIXED 
A 2 ( N ) =P,P n 
A3 t 1 ) =0.D0 
A1(N)=1.PP 
A'l ( 1 ) = 1 • DO 
A3 ( N ) = o . DO 
A2 ( 1 ) = ^ » dp 
DP 99 J=1,N 
AA( 1 » J ) = p , 

99 A A ( N > J ) = o • 

C PFADl IS VALUE OF FUNCTION AT X=0 
Ft 1 ) = R F AD 1 
G(1)=F(1) 

C PFAD2 IS V A L U c OF FUNCTION AT X = L 
F ( N ) = R F A O 2 
C-( N I - F ( N ) 

I TFR = 0 

C BEGIN ITERATION CYCLF 
115 DP 1 30 J=] ,N 

SU k ' = G ( J ) 

DO 133 JJ=1,N 

133 CI|*’rSMK-JA( J, JJ)*F{ JJ) 

1 *3 O Y ( J ) = c 1 IM 

WWt l) = A^n)/Al(l) 

GG ( 1 ) = Y ( 1 I /A] ( 1 ' 

DO 135 J=2»N 

GGt J ) = ( Y ( J 1 - A 2 ( J ) *GP- (J-l) )/(Al (J) — A2(J) *WW( J-l ) ) 
1 3 c WW t J ) = Ai ( J ) / ( A 1 ( J ) -A2 ( J ' *WW ( J-l ) ) 

nnl=n 

SAV-F ( N ) 

F ( N ) = G G { N ) 

roc-DAPp ( ( g a v — r ( n ) ) / R ( n ) ) 
x = m-1 

no 136 J=2»N 
5 AV = F t K ) 

F ( <) =GG ( < ) -WW ( K ) *F ( X + 1 ) 

A-PABF ( ( SAV-F ( X ) ) /F ( X ) ) 



IF(A.LE.EPS) GO TO 136 
EPS = A 
NNL = K • 

136 K = K-1 

I TER= I TER + 1 
C CONVERGENCE TEST 

IFtEPS.LE.EP) GO TO 160 
151 IF( ITER.LF.IT )G0 TO 115 

160 V(101)=0. 

RB = 0 • 

DO 170 J= 1 ♦ 21 

RA=RB**2 

RC= 1 . -RA 

RD=1 .+RA 

FRA=-4.*RA 

DO 180 L= 1 *100 

XL=L-1 

XL= LBAR-XL*H 
XR=XL**2 

Y ( L ) =F ( L )*XL* ( XR+RC )/ ( ( XR+RD)**2+FRA)**1,5 

180 CONTINUE 

51 = 0. 

52 = 0 . 

DO 181 I =2 .N.2 
Sl = Sl + Y ( L ) 

181 S2=S2+Y(L+11 

WW( J ) =H3* ( Y-Y (101)+4.*S1+2.*S2)*RB 
170 RB = RB+ .05 

51 = 0. 

52 = 0 . 

DO 182 J= 2 » 20 » 2 
?l=Sl+WW( J ) 

182 S2=S2+WW( J+l ) 

D IND=.2#C*( WW-WW ( 2 1 )+4.*Sl+2.*S2) /3. 

DTH1=(F(2)-F)/H 
DTHL= ( F ( N ) -F ( N-l ) ) /H 
PG=PI ND+PD 
DO 164 J= 1 >27 
L= I NUhX ( J ) +1 
164 Y(J)=F(L> 

PS=2.*D*C*DTHL 
C PRINT ANSW C PS 

WRITE (6 .101 ) LBAR.D.C.RFADl .RFAD2 .0TH1 ,DTHL 
WRITE (6.102) ITEP.FPS, PIND.PG.PS, H 
WRITE (6,103 ) ( INDEX! J) ,Y(J)»J=1»27) 

Y ( 28 ) = DTH 1 
Y ( 29 ) =L> l HL 
Y(30)=ITFR 
Y ( 3 1 ) = EPS 
Y(^2)=PIND 

CALL BCDUMP ( Y ( 1 ) »Y(32) j 
WR I T E ( 6 . 104 ) 

WRITE (6, 104) 

GO TO 1 

101 FORMAT ( 6H LBAR = F7 . 4 , 3H ♦ D= F 1 5 . 8 » 3H , C = F 2 . 0 . 9H . THET A 0 = G15.8, 

* 9H, THETA L = G1 5 . 8 » 8H , DTHDX 0=G 1 5 . 8 > 8 H , DTHDXL = G1 5 . 8 ) 

102 FORMAT ( 6H ITER=I3»5H,EPS=F8.6»6H,PIND=G14.8»4H',PG=G14.8» 

* 4H.PS=G15.8»8H,H=F8.6 ) 

103 FORMAT! 1H09( 14. 1PE10.3J/1H 9 ( 1 4 , 1 PF 1 0 . 3 ) / 1 H 9 ( I 4 , 1 PE 1 0 . 3 ) ) 

104 FORMAT ( 1HK ) 

105 FORMAT (6F10. 5 ) 

EMD 
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TABLE II. - INDIRECT TRANSMISSION FACTORS 
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TABLE III 


MATCHED BOUNDARY CONDITION SOLUTIONS 
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Figure 1. - Parameters useful for solution of upstream and downstream wall expressions (eqs. (21) and (27)) 
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Figure 2. - Variation of equilibrium arrival rate with temperature for various metal vapors (ref. 20 
and eq. (41)). 
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Figure 3. - Variation of adsorption time with energy for various 
temperatures. Constant, t q = 10“ 13 second (eg. (C7)). 
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Figure 4. - Variation of direct transmission factor with tube 






"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the. atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof .” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless 

of importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distri- 

bution because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in con- 

nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 

and initially published in the form of journal articles. 

SPECIAL PUBLICATIONS: Information derived from or of value to 

NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may he obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



